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1 Introduction 



The quantum Hall effect (QHE) describes the generic behavior of gapped fermions 
in 2+1 dimensions with broken parity]^ A typical physical manifestation involves 
a low-temperature electron gas confined to the interface of a heterojunction in a 
strong, transverse magnetic field. As a result of its striking and robust experimental 
signatures, the QHE has been a vigorous area of research since its discovery three 
decades ago. 

Much of the recent focus has been on the fractional quantum Hall effect (FQHE), 
in which the filling fraction z/, the ratio of the charge density to the magnetic field 
in units of e/h, takes certain fractional values. One reason a full explanation of the 
FQHE remains elusive is that it results from strongly-coupled dynamics, as evidenced 
by interesting features such as fractionally-charged quasiparticles (2||3]. 

The integer quantum Hall effect (IQHE), where i/ G Z, was discovered first because 
it is much easier to detect experimentally. It is also much better understood since, 
for the most part, it can be understood in terms of free quasiparticles organized into 
Landau levels. However, recent interferometry experiments ji] indicate that perhaps 
the situation is not so simple and interactions are indeed important. 

Gauge/gravity duality has proven to be a powerful tool for investigating strongly- 
coupled systems, such as those arising in many interesting condensed matter systems. 
The QHE, in particular, has proven to be amenable to holographic study. For ex- 
ample, phenomeno logical "bottom-up" models [5]-|9] of dyonic black AdS black holes 
have been proposed to capture certain observed features of QHE physics, in particular 
those related to SL{2, Z) duality, but do not include a mass gap. 

Brane constructions realizing quantum Hall physics 10 12 predate the current 
holographic approach. In more recent brane models 13 21 , 2+1-dimensional brane 
intersections with ij^ND = 6 yield a low-energy spectrum of fundamental fermions]^ 
These models are related to the Sakai-Sugimoto model of QCD in 3+1 dimensions [23] . 

The philosophy of this type of "top-down" approach is essentially pragmatic. 
Starting with a tractable string theory solution yields a holographic dual of the known 
field theory whose properties can be studied and which can hopefully then be matched 
to interesting physical systems. Previous models have been based in different ways 
on D3-D7 intersections, but there are other possibilities which should be analyzed 
in order to determine which features are robust and universal across models and the 
degree to which the others vary. 

With this approach in mind, we construct a holographic model of N species of 
2+1-dimensional fermions coupled to a strongly-interacting 2+1-dimensional gauge 
field using a probe D8-brane nontrivially embedded in a thermal D2-brane back- 
ground. We refer to this model as the D2-D8' model with the prime denoting the 



^For a review see, for example, l|. 

^See also 22 for a M5-brane model of the 1+1-dimensional edge of a QH system. 
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nonstandard embedding and emphasizing the analogy with the D3-D7' model of [20] . 
The embedding is stabilized by a worldvolume flux wrapping an internal cycle. The 
charge density and transverse magnetic field are holographically encoded in additional 
components of the worldvolume gauge field. 

We find that the two classes of embeddings correspond to different phases. Min- 
kowski embeddings (MN), for which the D8-brane smoothly caps off before entering 
the black hole horizon, represent gapped QH states. The dynamically-induced mass 
gap for charged excitations is roughly given by the distance between the tip of the 
D8-brane and the horizon. Black hole embeddings (BH), where the D8-brane crosses 
the horizon, are ungapped, metallic states. For a given charge density, there is a 
single MN embedding; this solution has filling fraction per fermion species u = 1, 
meaning that the D2-D8' system is a holographic model of the IQHE. 

We compute the electrical conductivity and the phase diagram and find some 
qualitative agreement with experimental observations. As expected in the QH state, 
the Hall conductivity a^y = ^ and the longitudinal conductivity a^x = 0. The QH 
state exists only at low temperature; we find a first-order phase transition rather 
than a crossover to a metallic state when the temperature is as large as the mass gap. 
Moving away from the QH state by changing the charge density or the magnetic field 
again leads to a phase transition to the metallic state rather than a crossover. The 
conductivity of the metallic state qualitatively obeys the semicircle law [8 24 26], 
albeit in a very rough sense, tracing out something like a hyperbola in the {axy, Cxx)- 
plane as the magnetic field is varied. 

Some discrepancies with real QH systems are due to unphysical simplifications of 
the model. In particular, the lack of impurities, which are required for the charac- 
teristic plateaux in the Hall conductivity, is the reason that the MN embedding only 
exists, for a given charge density, at a single value of the magnetic field. In addition, 
the conductivity of the perfectly clean metallic phase has unphysical behavior as zero 
temperature is approached. 

In many respects, this model reproduces many of the features found in the D3-D7' 



model of 20 . One important difference is that while the D2-D8' model is dual to a 
fully 2 -fl- dimensional boundary theory, the D3-D7' describes a defect theory in which 
fermions are confined to 2+1 dimensions while the gauge bosons propagate in 3+1 
dimensions. Another significant difference is that rather than exhibiting the IQHE, 
the D3-D7' models the FQHE with a filling fraction set by the internal fluxes. 

In the following section, we describe in detail the D2-D8' system, and in Section [3} 
we find the D8-brane embeddings which correspond to the quantum Hall and metallic 
states. In Section |4] we compute the conductivity of these two states, and in Section 
[5] we analyze the phase structure of the model. We conclude with a summary and 
discussion in Section [6l 
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2 The D2-D8' system 



In order to construct a nonsupersymmetric system of fermions in 2+1 dimensions, 
we consider a setup of N D2-branes and a single D8-brane with ij^ND = 6. Internal 
flux is needed to stabilize this nonsupersymmetric system. The low-energy spectrum 
of bifundamental strings in a jj^ND = 6 system contains only charged fermions and 
no charged bosons. In this case, the boundary field theory is 2+1-dimensional SYM 
coupled to N species of charged fermions. Similar constructions have been used, for 



example, in the D3-D7 15 and D3-D7' 20 models 



We will work in the limit iV ^ 1, so we can treat the D8-brane as probes in the 
background sourced by the D2-branes. We will also take the usual 't Hooft limit, 
where QgN 1. 

2.1 D2-brane background 

We consider the near-horizon background of thermal D2-branes|^ The 10- dimensional 
metric is 

L-^dslQ = u^ {-f{u)dt^ + dx^ + dy'^) +u-U--^ + u'^dnle] , (1) 

where = Gn'^gsNl^ and the thermal factor f{u) = 1 — {^)^ ■ Note that this space 
is not asymptotically AdS. The corresponding Hawking temperature is T = ^j^u'^'^. 
Furthermore, the background dilaton and the RR 3-form potential are 
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gs (^^^ ' , C(^) = -u^L^dt AdxAdy . (2) 



We will work in dimensionless (lowercase) coordinates, which are related to the 
physical dimensionful (uppercase) coordinates as: x^ = X^/L, u = U/L. In general, 
lowercase Latin letters will denote dimensionless quantities, and uppercase letters will 
denote the corresponding physical quantities. 

We choose to parametrize the in a somewhat unconventional way, analogous 
to the D3-D5 system of p~5| and the D3-D7' model of 20 , as an 5^ x fibered over 



an interval. Introduce the angles ip, 6i, 0i, 6*2, 02, and ^ such that: 

dQla = dip'^ + sin^ ip [dOl + sin^ 6*1^01) + cos^ if) (rf^^ + sin^ S,del + sin^ ^ sin^ 02d(t)l) . 

(3) 

The angle if) ranges from to 7r/2, the angles ^, ^i, and 6*2 range from to vr, and the 
angles 0i and 02 range from to 27r. 



See 27 for a thorough discussion of the D2-brane background and the regimes of vaUdity of 



various descriptions. 
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We will need the RR 5-form potential C-^-* dual to C^^\ Taking the Poincare dual 
of the 4-form field strength G^'^^ = dC^^^ gives 

G^^^ = sin^ ^p cos^ i/j dip A dfis^ A dils^ . (4) 

This can be integrated with respect to ip to find 

C(5) = c{tp)L^dQs^ A dQs^ , (5) 

where c{tp) = | (sin?/' — | sin(3^/') — ^ sin(5^/')). We have fixed the constant of inte- 
gration by choosing the gauge c(0) = 0. 



2.2 D8-brane probe 

The probe D8-brane fills the spacetime directions t, x, y, and the radial direction u 
and wraps the 5*^ and the S^. The embedding in the -^-direction will be a function 
of u, as ip{u). The induced metric reads 

L~'^ds\)^ = {-fdf + dx'^ + dy"^) + M"i Q + m^^'^^ du^ 

+u~^ sin^ iidVtl + -u^^ cos^ iidVtl . (6) 

We would like to consider a boundary system at nonzero charge density in a 
background magnetic field. Although the boundary U(l) field is non-dynamical and 
the corresponding current is global, we can still mimic the effects of a background 
field on a local charge current]^ Holographically, this is accomplished by turning on 
components of the worldvolume gauge field; specifically, a constant field strength in 
the spatial directions corresponds to a transverse magnetic field, and nonzero charge 
density is dual to a radially varying electric fieldj^ 

2TTa'F^y = h (7) 
27r«'F„o = a'^{u) . (8) 

In addition, as we will see, in order to stabilize the embedding, we will also turn on 
a magnetic field wrapping the internal S"^: 

2TTaF = hL^ dVL2 . (9) 

Turning on this wrapped fiux can be alternatively viewed as dissolving smeared D6- 
branes into the D8-brane. 



^There is an extensive literature on this issue; for example, see 
^We choose a gauge where a„ — 0. 
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The DBI action for the D8-brane is then 

Sdbi = -/is j d\ e-'^^- det{g^, + 2na'F^,) (10) 

= -Af j duui coii^ip^il + /m27/>'2) [h^u + sin^) + ^ , (H) 

where J\f = Svr^Tgf and f 3 is the dimensionless volume of spacetime. As a resuh 
of the charge density and the magnetic field, the Chern-Simons action has a nonzero 
component: 

= -^{2TTa'f j C^^^AFAF = N' J duc{^)ha'o . (12) 
In order to make the Chern-Simons action invariant under gauge transformations of 



C^^\ we also need to add the following boundary term 20,31 , which affects the 
energy but not the equations of motion: 

SMry = -y (2W)2 J d (C(^) A F A A) = -Mc{i,)hao\Z.^ . (13) 

Since we have chosen the gauge c{il) = 0) = 0, there is no contribution at m = 00 and 

so Sbdry = ■N'c{Umin)hao(Umin)- 

From this action, we compute the equation of motion for ip: 

du I u^giu) I 1 H — ? I ^' I = —TT cos^ -ip sin -ip (2 cos^ sin^ ip — 3 (b'^u + sin^ -ip) ) 
V V ^ / / 9{u) 

-a^c(^)K , (14) 



where we have defined 



g{u) = I = . (15) 

;i + - (1 + S) 



1 + — ) a'o = - hc{ij) = d{u) , (16) 



The action is independent of oq, so the equation of motion for can be integrated 
once to give 

»(i + ^i" 

where d is the constant of integration and d{u) is the radial displacement field. In 
terms of the boundary theory, d is the total charge density, while d{u) is the charge 
due to sources located in the bulk at radial positions below u. 

In order to decouple the equations and express the ip equation in terms of the 
charge, we can solve for a'g in terms of d, obtaining 



°° ^' d■' + u^cos^^P{b''u + sm^^P){l + ^) 
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We can now write (15) as 



/ /J2 + m5cos6^/> (1 + ^) (62M + sinV) 
^ = rT^V 1 + fu^ij'^ ' ^^^^ 



and (14) can be rewritten as 



du\u^g\l-\ — ^ j ip' j = cos^ sinT/) (2 cos^ sin^ t/) — 3 + sin"^ ^/i)) 



—5 cos^ -0 sin^ '0 . (19) 

91 + ^ 



2.3 UV asymptotics 



Because supersymmetry is completely broken in this system, one might expect that 
the D8-brane embedding would be unstable. In flat space, a D2-brane and a D8- 
brane which are parallel with jj^ND = 6 repel each other, leading one to anticipate 



a similar instability after taking the near-horizon limit 13 . Rather than analyzing 
the fluctuation spectrum directly, we will take the simpler approach of computing the 
dimension of the operator dual to ip to see if it has a nonzero imaginary part. The 
details of the calculations are presented in Appendix |Aj 

As expected, we see a tachyonic mode in the D8-brane embedding when there is 
no flux wrapping the internal S'^; the D8-brane slips to one of the poles of the S^. 



However, as in other similar systems, such as the D3-D7 model of 15 and D3-D7 



model of 20 , worldvolume fluxes can stabilize the probe brane. For any nonzero b, 
there are apparently stable D8-brane embeddings. Of course, a detailed analysis of the 
fluctuations will be required to guarantee stabihty, but at least the most dangerous 
tachyon has been removed. 

This result is in contrast to what has been seen in the analogous D3-D7 systems 
15 ,20 , where a minimum nonzero flux was required]^ Also unlike the D3-D7' system, 
both the asymptotic angle, V'oo = 0, and the subleading behavior, A = — 1 or —3, are 
fixed and are independent of the size of the flux. 

In the flat space embedding, the fermion mass is given by the minimal distance 
between the D2-brane and the D8-brane. As in the D3-D7' system, in the near-horizon 
limit, we identify the coefficient of the leading term in the asymptotic expansion of 
ip with the mass and the coefficient of the subleading term as the chiral condensate: 

^ = - - 4 + ^(«"') • (20) 



Note, however, that in the D3-D7' system of 20 , for the special cases of -000 — and tt/2, there 



is no lower bound on the flux required for stability. 



7 



3 Embeddings 



Having understood the UV limit of the D8-brane embedding, we now consider the 
full embedding. There are two classes with different IR behavior, black hole (BH) 
embeddings, which enter the horizon at mt, and Minkowski (MN) embeddings, which 



smoothly cap off at some Uq > ut- As argued in 20 , the MN embeddings holographi- 
cally reproduce the properties of a quantum Hall fluid, while the BH solutions exhibit 
the metallic behavior characteristic of the regions between quantum Hall plateaux. 

For the figures presented here, we have chosen values of the parameters ut, m, 
d, h, and b for illustrative purposes. The scaling u — )■ Am, ut Xut, rn — )■ Am, 
d X^/^d, h X^/^h, and b X-^/% is a symmetry of the equations of motion, and 
so there is, in fact, only a four-dimensional parameter space. We find that qualitative 
results are independent of these choices. However, we can not conclusively rule out 
the possibility of different behavior in particular corners of parameter space. 

3.1 MN embeddings 

For the D8-brane to smoothly end at Mq, one of the internal spheres must shrink 
to zero at the tip. Because of the flux wrapped on the S*^, we choose the to 
shrink instead, giving ?/'(mo) = vr/2 and as u — > mq, ip' — )■ — oo. As in the D3-D7' 



model 20 , because the tip of the D8-brane cannot support a worldvolume charge, all 



the charge must be induced by the magnetic field via the Chern-Simons term (12). 
This implies that for MN embeddings, the ratio of the holographic charge density d 
and the magnetic field h is fixed: 

^=c(./2) = !. (21) 

The filling fraction per fermion species v is defined by the ratio of the physical charge 
density D = ^^^^^jj^d to the physical magnetic field H = divided by the number 
of species A^. We find that 

2nD 3d 

V = = — = 1. (22) 

The MN embedding is therefore a holographic model of an integer, rather than frac- 
tional, quantum Hall fluid. 

Alternately, we can compute u from the 3-dimensional effective Chern-Simons 
term 

Scs = -^T8(2W)2 j ^(6) AAAF . (23) 

Because our MN embedding goes from ip = to n /2, we integrate over the whole 5*^ 
to find 

= 87r^L^ / d^c'i^) = {2nhfgsN . (24) 
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Plugging this in to (23) gives: 



cs 



N 

4:71 



AAF . 



(25) 



The level k can be read off from the coefficient of the Chern-Simons term f-; we see 



A- 

47r' 



directly that k = N. And since the level k is also given by the total filling fraction 
Nu, this gives the same result as above. 
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Figure 1: The mass m as a function of the position of the D8-brane tip uq, for 
b = 1, d = 1, and h = 1.5. The bottom, blue curve shows ut = 0.1, the middle, 
black Ut = 0.666, and the top, red ut = 1-2. Notice that the minimum mass mmin 
increases with temperature. For m = 0.1, at the low temperature there is one MN 
solution, two solutions at the intermediate temperature, and then no solutions at the 
high temperature. 



To find the complete embedding, we solve equation (19) numerically by shooting 
out from the IR tip and extracting m and c from the UV asymptotics. 

For certain initial conditions, ipi^^) becomes negative at some u and then asymp- 
totically approaches ip = from below, resulting in m < 0. While there is nothing a 
priori wrong with a negative fermion mass, in this case it is a sign that the low-energy 
approximation is breaking down. When, at some finite u, ip{u) goes to zero with fi- 
nite the D8-brane self-intersects at a conical singularity. For such a configuration, 
massive modes of the D8-brane, which are excluded from the DBI action, become 
light and, in fact, tachyonic, and the embedding is unstable. As a result, we only 
consider positive masses. 

We show the dependence of m on uq for several temperatures in Fig. [TJ For a 
given temperature, and fixed b and d, there is a minimum mass m^m for which a 
solution exists. At temperatures sufficiently below the mass gap, i.e., ut <^ Uq, this 
minimum mass is negative, and there exists a unique solution for any given positive 
m. However, nimin increases with the temperature, becoming positive, and eventually. 
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when m^m ~ 'm, a second MN embedding at smaller uq appears. As we shall see 
in Section |5| this embedding is thermodynamically unstable. We also expect it to 
be perturbatively unstable as well. As the temperature is increased further, the two 
solutions merge and disappear; at sufficiently high temperature, there are no MN 
embeddings with mass m because rrimin > 
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Figure 2: The mass gap in the MN embedding, mg^p = Uq, as a function of the 
magnetic field h for ut = 0.01, 6=1, and m = 0.1. The blue curve is a fit to 



A characteristic property of a quantum Hall fluid is a mass gap, and, as shown 



m 



20,32 , MN embeddings holographically realize the mass gap, denoted by m. 



go.pi 



by avoiding the horizon. Although we postpone a detailed analysis of the fluctua- 
tion spectrum for a future work, we can easily estimate the mass gap for charged 
excitations. Charged quasiparticles are holographically dual to strings stretched from 
the horizon to the D8-brane tip, and at zero temperature, the mass of such strings, 
and thus the mass gap for charged excitations, is simply ntgap = Uq. For an integer 
quantum Hall fluid, the mass gap increases linearly with the magnetic field; here, 
however, the gap scales as h^'^^ to very good approximation, as shown in Fig. M 



3.2 BH embeddings 

Black hole embeddings are those which enter the horizon. As with MN embeddings. 



we solve (19) numerically by shooting, in this case from the horizon. We need only 



to choose a single IR boundary condition, the location where the D8-brane intersects 



the horizon i/jt, because equation (19) becomes ffist-order in the near-horizon limit. 



Some initial conditions lead to embeddings with negative mass, which, as discussed 



in Section 3.1, we disregard as unstable artifacts. 
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Figure 3: The mass m as a function of iIjt-, with h = 1, d = 1, and h = 1.25. Curves 
for three temperatures are plotted: bottom, blue for = 0.1, middle, black for ut = 
0.53, and top, red for ut = 1-2. For m = 0.1, for example, at the low temperature 
there is a single solution at ipT ~ 1.44. There are three solutions ipT ~ 0.25, 0.88, 
and 1.10 at the intermediate temperature, and at the high temperature, one solution 
at ~ 0.07. 

For generic h and d, the behavior of m as a function of ipT is as shown in Fig. [3j At 
low temperature, there is a single solution at ipT ^ 7r/2 for any given positive mass. 
But, as the temperature is raised, two embeddings appear; as we will see in Section 
[5| the one at smaller ipT is stable, while the other is not. At still higher temperature, 

n/2—tffT 

usini/f 1.00 f 




Figure 4: In the left panel, as the magnetic field approaches locking, h ^ ^d, the 
large ipT BH embedding (dashed pink line) becomes increasingly spiky, ipT ~ 7r/2, and 
approaches the MN embedding (solid blue line). Both embeddings are for m = 0.1, 
b = d = 1, and ut = 0.3. The MN embedding has h = 1.5, while the BH embedding 
is /i = 1.49. The black hole is indicated by the quarter circle near the origin. In the 
right panel, the approach to a MN embedding, also for m = 0.1,b = d = l,is shown 
as a log-log plot of tt/2 — ^/jt versus |A/i|, illustrating | — ^/'t ~ \Ah\^^'^. The solid 
line is for Ah > 0, and the dashed line for Ah < 0. 
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the unstable solution merges with the large ipT solution, and they disappear, leaving 
only the small ipT embedding. 

As in [20], as the magnetic field is tuned such that z/ — )■ 1, the large ipT BH solution 
starts to resemble a MN embedding connected to the horizon by a thin spike. Fig. |4] 
shows how — )■ 7r/2 as the magnetic field is tuned to u = 1. We find that for small 
Ah = h- |d, 

I _ |A/i|i/2 . (26) 

Although the transition from a BH to a MN embedding at Ah = appears 
discontinuous, certain physical quantities, such as the conductivities are continuous, 
as we will see in Section |4} We will investigate more closely in Section 5^ and find a 



phase structure analogous to the MN to BH transition in the supersymmetric D3-D7 
system 



33,34 



4 Conductivity 

To compute the conductivities in both the BH and MN states, we follow the well- 



known Karch-O'Bannon procedure 35 36 . We holographically introduce an electric 



field in the x-direction and currents in the x- and y-directions, as follows: 

27ra'Fo. = e (27) 

27ra'Fu, = < (28) 

27Ta'Fuy = a'y . (29) 

The DBI action is now 

Sjjsi = j d-uu^ cos^i^^Jv (62M + sin^V^) , (30) 



where 



5 i ^'^ 



+K' + /^l-7^J<-2«o<^- (31) 
There is also an additional term in the Chern-Simons action 

Scs = U I du c{ip) [ha^ + ea'y) . (32) 



The equations of motion for the gauge fields can be computed and integrated 
once, introducing constants of integration d, jx, and jy as follows: 

7U^ + ^J^° + ^" 



d^d-hc{^) = '-,{{l + -E)ao + -^a'y] (33) 



jx = ga'x (34) 



12 



where the function g is now 



g = fu'/^ cos' ^ . (36) 

Using some algebra, we can solve for the gauge fields in terms of the currents: 



/m5 / u^'yj V X 



a. 



IH-^li-S^-^yii^^. (39) 



where 



X = f (^1 + ^^- ^^(^fu' cos'' ^{b'u + sin' ^P)+fd'-jl -J 
fhd + ejy 



(40) 



4.1 Quantum Hall fluid 

For the MN embeddings, we use a modified version of the Karch-O'Bannon method 



20 ; the conductivities can be computed from the requirement that the gauge fields 



be regular at the tip mq, which is accomplished by 

d = 3. = ~3y = ^ ■ (41) 

The vanishing of d{uQ) implies 

^ = c(vr/2) = 2/3 . (42) 

Because = 0, there is no longitudinal current and 

a,, = . (43) 

Lastly, the vanishing of jy gives jy = ec(7r/2), from which we find the physical Hall 
conductivity 

_Jy _ (27raOW j, _ AT 

- E L^Vs e 27r ' ^ ' 

Previously, in Section 3.1 we found the filling fraction per fermion species u = 1. 
Therefore, we find 

Nu 
2^ 
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which is the expected transverse conductivity in a quantum Hall state with fermions, 
each with filling fraction v. 

It should be noted that the conductivities computed here are independent of the 
temperature. The tip of the D8-brane is away from the horizon, and therefore the 
condition that the charge at the tip vanishes should not depend on the horizon's 
exact position. In particular, in terms of the two-component description of the current 
carriers 20 37 , all of the charge is of the dissipationless type, smeared along the radial 
direction of the D8-brane. The strings stretching from the horizon to the D8-brane 
are massive, and fiuctuations producing charge-anticharge pairs are gapped]^ 

We have found a vanishing longitudinal conductivity. However, this method only 
takes into account massless charge carriers. The QH state is gapped, and the lightest 
charged excitations which could support a current are thermally pair-produced strings 
stretching from the horizon to the D8-brane. Because such strings have masses of the 
order of mgap, the Boltzmann factor for such fiuctuations is of the form e~'-^™'9''p/-^, 
for some unknown constant C . Including such suppressed thermal pair production 
would likely yield a^x ~ eT^'^^^^l'^ . Such an exponentially-suppressed conductivity is 



observed in real quantum Hall systems 38 , 39 



4.2 Metal 

For the BH embeddings, we can compute the longitudinal and Hall conductivities 



using the standard Karch-O'Bannon method 35 ,36 ; the currents are fixed by the 
condition that the action be real. In particular, this requires that X have double 
roots at the pseudohorizon m^,, which implies the following three conditions: 

= 1 + ^-^ 46 

ejy* = f*hd^ (47) 
Jl+Ty. = f.uUos^^p4b\, + sm^^P,)+fJl. (48) 

While these conditions, in principle, give the full nonlinear conductivity, we will 
specialize to the case of small electric field to obtain analytic expressions for the linear 



response. The first condition (46) gives the location of the pseudohorizon which, to 
leading order in e, is 



e2 



The thermal factor at the pseudohorizon, to leading order in e, is then 



^This somewhat intuitive explanation of the mass gap has been verified in the D3-D7' model by 



a more careful fluctuation analysis 32 . We will leave the study of the fluctuations of the D2-D8' 
system for a future work. 
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The second condition (47) gives the physical, hnear Hall conductivity 



a 



xy 



Air 



d + c{'lljq 



(51) 



and the third condition (48), when combined with the condition for the Hall current 



(47), gives the physical, linear longitudinal conductivity 



(Tr. 



u 



5/2 
T 



d? + cos6 {bW + sin^ ^jt) + h'^) ■ (52) 



We numerically evaluate the conductivities and display the results in Fig. [5j Both 
the longitudinal and Hall conductivities are nonzero, as one would expect for a metal- 
lic phase. Note that as d/h 2/3, the conductivities smoothly approach their QH 



values, (43) and (45). As was seen in both the D3-D7' model pO] and the Sakai- 



Sugimoto model 37 , the Hall current is made of contributions from two types of 



current carriers; one with charge density d is located at the horizon and loses energy 
to the heat bath, while the rest of the charge density, hc{ipT) = d — d, is distributed 
along the D8-brane above the horizon and suffers no dissipation. 




0.00 
0. 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 



Hj- 



Figure 5: For the metaUic phase, the conductivity is plotted as a function of temper- 
ature with h = 2, m = 0.1, and b = d = 1 in units of the Hall conductivity in the 
QH state (J^y^- On the left, the longitudinal conductivity a^x scales as u^'^ at low 
temperatures. On the right, the Hall conductivity a^y is approximately independent 
of temperature. 



One seemingly non-metallic property is the low-temperature behavior of the longi- 
tudinal conductivity. From (52) we see that, for vanishing magnetic field, a^x diverges 

—5 /2 ^ 

as Urp . This is due to the unbroken translational symmetry of the D2-D8' system; we 
have introduced neither impurities nor edges, both of which would break translational 
invariance. Similarly, (52) implies that when h is nonzero, a^x vanishes as u^"^, which 

the full Lorentz invariance is restored and 



is a consequence of the fact that at T 
the longitudinal conductivity must vanish. For a real metal, a^^ approaches a finite. 
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. QH 




0.000 
0.0 



. QH 



Figure 6: The longitudinal conductivity is plotted against the Hall conductivity 
a^y for varying magnetic field h in units of the Hall conductivity of the QH state 



aib = d = 1 and m = 0.1. The temperatures are, from top to bottom, ut = 0.5 (solid 
black), 0.4 (dashed blue), 0.3 (dotted red), and 0.2 (dot-dashed pink). For h = 1.5, 
at the MN embedding, [o^yj a^^ .a^^j a^^) = (1,0). The curve proceeds to the left 
as h increases, and for h — )■ oo, the curve approaches the origin. 



nonzero value at small temperatureqj due to scattering off impurities. Here, charge 
carriers are only scattered by thermal fluctuations and so behave at zero temperature 
as free particles. 

In a quantum Hall system at low temperature, as the magnetic fleld is varied from 
one plateau to another, the conductivities sweep out a semicircle when plotted on the 
(axy, axx)-p^dJie 24 . The D2-D8' model contains only a single QH state and therefore 



strictly speaking can not model the transition between such states or reproduce this 
semicircle behavior. However, in a system exhibiting only the IQHE, as the magnetic 
fleld is increased beyond the u = 1 state, the system becomes a quantum Hall insulator 
and the conductivities trace a semicircle which approaches the origin, axx = (^xy = 0, 



at large magnetic fleld 25 26 . As illustrated in Fig. 6 the D2-D8' model qualitatively 



reproduces this behavior. One important discrepancy, apparent in (52), is that, rather 
than reaching a maximum at axx = 1/2, the maximum of the "semicircle" scales at 
low temperature as u^'^. We interpret this as a probable consequence of the unbroken 
translational invariance. 



5 Phase structure 

Our flrst task in order to investigate the phase structure of the D2-D8' system is 
to compute the free energy. The regulated, on-shell Euclidean bulk action gives the 

^This assumes, of course, that there is no transition to a superconducting phase. 
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grand potential of the boundary theory: 

T, h) = j^S^_gi^^ii = jj: {SdbI+CS + ^bdry + Sct) 



shell ' 



where is the chemical potential and S^ry is the boundary term given by (13). We 
regulate the divergences of the action by adding a gauge and diffeomorphism invariant 
counterterm0 

ScT = - ^^""f j dT. e-^y/det(7^, + 27r«%,) (53) 

which is integrated over the UV cutoff surface S with induced metric 7. By choosing 
the cutoff to be the surface u = uuvi the counterterm is then 

S^^ = -U\h\'^{l + 2i,ly) (54) 

which cancels the divergences but also contributes a finite piece proportional to |6|mc. 

Since we are working in the canonical ensemble, we need to perform a Legendre 
transform to obtain the free energy: 

F{d,T,h) = n{fi{d),T,h) + dfx . (55) 

We can write the chemical potential fi(d) = ao(oo) = dua'n + ao(umin), such that 



F = j du cos^i)^{l + /m2^'2 - a'^) [h'^u + sin^ ^) + ^ 

+ rfa'o - c{i))ha!^ + jj-Sct + dao{Umin) - c{i){Umin))hao{Umin) ■ (56) 

For BH embeddings, ao(Mr) = 0, so the last two terms vanish. For MN embeddings, 
recall that d = d — c{7v/2)h, and so the last two terms cancel. The free energy is then 



F = J du cos^ (1 + fu^ij'^ - (a'o)2) (b^u + sin^ ^) (^1 + + da'^ 

+^ScT . (57) 

We can now compute the free energy for the various embeddings and, by comparing 
them, determine the thermodynamically preferred state as a function of the external 
parameters. 



''Any such counterterm with these symmetries and which cancels the divergences will suffice. 
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5.1 High-temperature phase transition 

We begin with the special case where the magnetic field and charge are locked, 
d/h = 2/3 and find a first-order phase transition from a MN embedding to a BH 
embedding when the temperature is increased. This critical temperature is approxi- 
mately the mass gap, m^** < uq. The free energy as a function of temperature, shown 
in Fig. [7| takes the standard swallowtail form; near the phase transition, there exists 
an additional unstable phase which, depending on the temperature, is either a MN 
embedding at uq close to ut or a BH embedding with ipT close to 7i/2. 
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0.75 



0.80 
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Figure 7: A plot of the free energy F as a function of the temperature ut in the 
locked phase, d/h = 2/3, for m = 0.1 and b = d = 1. The solid blue curve shows the 
stable MN embedding, the dashed black curve is the stable, small ipT BH embedding, 
while unstable MN solution is solid green and the unstable BH embedding is dashed 
red. The swallowtail structure indicates a first-order phase transition at ut ~ 0.6. 



A similar first-order phase transition was seen in the D3-D7' model 20 . In real 



quantum Hall fluids, however, when the temperature is of the order of the mass gap. 



the systems smoothly crosses over to metallic behavior 38 , 39 



For generic charge and magnetic field, the low temperature phase is no longer a 
MN embedding but instead a large ipT BH embedding. There is still a first-order 
phase transition as the temperature is increased, but now it is from a BH embedding 
with large ipT to one with small ipx- The phase diagram for the temperature-magnetic 
field plane is shown in Fig. [8j For fixed d, the critical temperature grows with h. The 
line of first-order phase transition ends at a critical point; for sufficiently large h, there 
is always only a single BH embedding, and ipr decreases smoothly with increasing 
temperature. 
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0.8 r 




Figure 8: Phase diagram in the temperature-magnetic field plane for b = d = 1 
and m = 0.1. The solid blue curve is a line of first order-phase transitions from 
large ipT spiky BH embeddings at low temperature to small ipT embeddings at high- 
temperature. This line ends in a second-order critical point {h,UT) ~ (2.3,0.69) 
indicated by a black dot. The dashed black vertical line indicates the locked phase, 
d/h = 2/3 where there is a MN solution rather than a large ipx BH embedding. 



5.2 Transition from the QH state 

There is also a transition from the QH state to a metallic state, that is from a MN 
solution to the large ipT BH embedding, as either the charge or magnetic field is 
varied away from d/h = 2/3. A similar transition from MN to BH phases occurs 
in the supersymmetric D3-D7 model, where the MN embedding only exists at zero 



density 33,34 



The free energy has a kink as the system passes through the MN solution; at fixed 

dF 
dd 



the magnetic field, the chemical potential /i = appears to be discontinuous across 



d = T^h, as shown in Fig. |9| The jump in /x can be calculated by noticing that for very 
spiky BH solutions, Oq = ±1 along the length of the spike, with the sign depending 
on whether d/h ^ 2/3 from above or below. Away from the spike, the limit of Oq 
as d/h — 2/3 is otherwise smooth. Since the length of the spike is uo — ut, we find 
AyU = 2{uq — Ut)- 

This discontinuity is an artifact of working in the canonical ensemble. The problem 
is that the MN solution exists at a single d but for any fi. Because the D7-brane avoids 
the horizon, oq does not have a Dirichlet IR boundary condition, and ao(c)o) = /i is 
unfixed by the dynamics and can take any value. As a result, the canonical ensemble 
is not well defined. 

We will switch instead to the grand canonical ensemble, where the nature of the 
transition can be more clearly understood. MN solutions exist for all /i, but there is a 
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d 
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Figure 9: The charge density d as a function of the chemical potential /i for fixed 
h = 3/2, b = 1, ut = 0.1, and m = 0.1. Away from d = 1, the embedding is a large 
ipT BH (blue curve), while at (i = 1 it is a MN embedding with mass gap uq ~ 0.83 
(dashed black line). The range in fi without a BH solution is of size Afi ^ 1.46. The 
plot on the right zooms in to show that for fi between 0.048 and 0.074, there are two 
BH solutions in addition to the MN solution. The dotted vertical line illustrates the 
Maxwell construction, by which we can determine that the first-order phase transition 
occurs at /i ~ 0.068. 



range in /j, of size A/x where there are no BH embeddings, In this particular range, 
the system can only be in the MN phase. But, where the two solutions coexist, the 
BH phase has lower grand potential fl and is therefore preferred. At the large-/i edge 
of this range, d is linearly increasing with fi and so is continuous as the system goes 
from the MN solution to the BH solution; the phase transition is therefore marginally 
second-order. However, at the small-/x end, as can be seen in the zoomed in plot 
in Fig. [9} there are two BH solutions for a narrow range of fi. The branch which 
connects to the MN solution has fj; < 0, meaning it is unstablej^ As /i is lowered, 
the system jumps from the MN solution to the stable BH solution at lower d, thereby 
undergoing a first-order phase transition. The critical /i at which the transition occurs 
can be found via a Maxwell construction, as illustrated in Fig. |9j 

There is a surprising asymmetry between the transition away from the QH phase 
toward larger d and to smaller d. In the supersymmetric D3-D7 model 33 34 , where 
the analogous transition occurs at d = 0, the transition could be either first- or 
second-order depending on the temperature. However, at a given temperature the 
transition to positive and negative d was always identical due to charge conjugation 
symmetry. Here, this symmetry has been broken by our choice of h, although it 
remains unclear precisely why these transitions should be of different orders. 



^°Note that A/i = 2{uq — ut) is equal to the mass of a quark plus the mass of an antiquark. 
^"'^In an upcoming work on the fluctuation spectrum of this model, we will investigate the nature 
of this instability. 
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6 Discussion 



We have presented a holographic model of a 2+1- dimensional, strongly-coupled fer- 
mion system at nonzero temperature, magnetic field, and charge density which ex- 
hibits quantum Hall behavior with integer filling fraction. The gapped QH states of 
this D2-D8' system correspond to MN embeddings of the D8-brane, which do not 



enter the horizon. We found in Section |3.1| that these solutions exist only at low 
temperature, ut < mgap, and for a fixed ratio of the charge density to the magnetic 
field, d/h = 2/3, which implies u = 1. The mass gap scales with the magnetic field 
as /i"-^^. The Hall conductivity of this state, computed in Section 4T, takes the QH 
form axy = while the longitudinal conductivity vanishes. 

At higher temperatures and at generic d/h, the D8-brane is in a BH embedding 
which corresponds to a metallic phase, similar to that found in the transition between 



two QH plateaux. In Section |4.2[ the conductivity was found to somewhat reproduce 
the semicirlce behavior seen in QH transitions. However, the D2-D8' system shows 
unphysical behavior at the low temperature due to unbroken translation invariance. 



In Section |5.1[ we found a first-order phase transition from the QH state to the 
metallic state with a critical temperature of the order of the mass gap. This phase 
transition persists away from the MN solution, for generic charges and magnetic fields, 
as a transition between two different BH embeddings. At large enough magnetic field, 
the line of first-order phase transitions ends in a critical point. 

When the charge or magnetic field is varied away from the MN solution, there 



is a phase transition, as seen in Section 5.2 In the grand canonical ensemble, the 



transition away from the QH phase is first-order when going to smaller /i and second- 
order when going to larger /i. 

Having constructed a QH state with z/ = 1, one might inquire whether this model 



might include other integer filling fractions as well. As we discussed in Section |3.1 
for a MN embedding all the charge must be induced via the Chern-Simons term, 
which is in turn directly related to how much of the 5*^ x 5*^ is wrapped by the D8- 
brane. This quantity is topological in that it only depends on the endpoints of the 
D8-brane embedding; since the tip must be at ip = and the UV asymptotics 
demand ipoo = 0, the entire internal manifold is wrapped and u = if^ 

Hypothetically, one could obtain higher-integer filling fractions from multiply- 
wrapped embeddings. Such embeddings, like the negative mass embeddings discussed 
in Section[3} would self-intersect and suffer conical singularities when ip passes through 
zero. At these points, the low-energy description of the D8-brane would break down 
as a result of tachyons, and the embedding would most likely be unstable. 

Acknowledgments 

^^This is in contrast to the D3-D7' model where ipoo can take a range of values, leading to filling 
fractions between approximately 0.7 and 0.8. 
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A UV asymptotic analysis 

We investigate asymptotic solutions for the D8-brane embedding ijj{u) at large u. 
Considering first the case without internal flux b = and taking the ansatz ip ~ 
ipoo + where we assume Re A < 0, to leading order O (^""^j; the equation of 
motion (19) is 

= — 3 cos^ ipoo sin^ V'oo + 2 cos^ ipoo sin tpoo (58) 
which implies three possible solutions for tpoo 

i/joo = 0, 7r/2, arctan I \/ - I • (59) 




We can immediately exclude the trivial constant solutions = or 7r/2, because 
those correspond to D8-branes with vanishing volume and energy. 
The equation of motion to order 0{u^^'^'^'^) gives 

— A(A + 7/2) sin^ 'i/'oo cos^ '?/'oo = 3sin^-?/'oo cos^'?/^oo(2tan-?/'oo — 3cot-?/^oo) 

— 2sin'?/'oo cos^'?/'oo(4tan'?/'oo — cot-i/^oo) -(60) 

For ipoo = 0, this equation cannot be satisfied, so this is not, in fact, a solution. 
At this order -ipoo = 7r/2 is still a solution, but it similarly fails at order (9(m^/^+^'^). 
Furthermore, it can be shown that any solution with tp^o = or 7r/2 must be constant 
and therefore degenerate. 



For the third solution ijjoo = arctan ^y^j, equation (60) reduces to 



A2 + ^A + 10 = , (61) 



which has the complex solutions 
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Thus, the 6 = embedding is unstable as expected. 

To stabihze the embedding, we now consider 6 > 0. The leading term is linear in 



b, meaning that ( 19 ) is asymptotically 

= — 36m'^ cos^ ^/'oo sin -^oo , (63) 

and the only solutions are ipoo = and n/2. 

First, let us look at the tp^o = case. To the next order, the equation of motion 



([19j) is then: 

A6A(4 + Ay+^ = -3A6M3+A + 2AW+3^ . (64) 

Because we have assumed ReA < 0, the second term on the RHS is subleading, so 
we find 

A^ + 4A + 3 = , (65) 

which implies A = — 1 or —3. Note that A is real, as demanded by stability, for all 
nonzero b. 

Now, we will consider the ipoo = 7r/2 solution. The asymptotic equation of motion 



(jigj) is then: 

- A%A{4: + 4Ay+'^ = -3AV+^'' + 2A''b-\'+'^ . (66) 

Again, the second term on the RHS is subleading, so we obtain 

- A^A{A + 4A) = 3u-^^ , (67) 

which has not solutions for ReA < 0; thus, the power law ansatz is not a valid 
solution for ipao = vr/2. In fact, it can be shown that any nontrivial embedding will 
fail to satisfy the equation of motion, and therefore, only the degenerate constant 
solution ip = 7i/2 has this asymptotic behavior. 
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